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Abstract 

< 

■ Once the critical temperature of a cosmological boson gas is less than the critical temperature, 

. a Bose-Einstein Condensation process can always take place during the cosmic history of the uni- 

q \ verse. Zero temperature condensed dark matter can be described as a non-relativistic, Newtonian 

H 

gravitational condensate, whose density and pressure are related by a barotropic equation of state, 
with barotropic index equal to one. In the present paper we analyze the effects of the finite dark 
matter temperature on the properties of the dark matter halos. We formulate the basic equations 



> 

| describing the finite temperature condensate, representing a generalized Gross-Pitaevskii equation 

00 ' 

£Sj ■ that takes into account the presence of the thermal cloud. The static condensate and thermal 

cloud in thermodynamic equilibrium is analyzed in detail, by using the Hartree-Fock-Bogoliubov 
and Thomas-Fermi approximations. The condensed dark matter and thermal cloud density and 
mass profiles at finite temperatures are explicitly obtained. Our results show that when the temper- 
ature of the condensate and of the thermal cloud are much smaller than the critical Bose-Einstein 
transition temperature, the zero temperature density and mass profiles give an excellent descrip- 
tion of the dark matter halos. However, finite temperature effects may play an important role in 
the early stages of the cosmological evolution of the dark matter condensates. 
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I. INTRODUCTION 



At very low temperatures, all particles in a dilute Bose gas condense to the same quantum 
ground state, forming a Bose- Einstein Condensate (BEC), i.e., a sharp peak over a broader 
distribution in both coordinates and momentum space. Particles become correlated with 
each other when their wavelengths overlap, that is, the thermal wavelength At is greater 
than the mean inter-particles distance I. This happens at a temperature T < 2irh 2 /mkBn 2 / 3 , 
where m is the mass of the particle in the condensate, n is the number density, and ks is 
Boltzmann's constant jl, 2]. A coherent state develops when the particle density is enough 
high, or the temperature is sufficiently low. Quantum degenerate gases have been created 

new lines of 



by a combination of laser and evaporative cooling techniques, opening severa 
research, at the border of atomic, statistical and condensed matter physics [l- 

The Bose-Einstein condensation process was first observed experimentally in 1995 in 
dilute alkali gases, such as vapors of rubidium and sodium, confined in a magnetic trap, and 
cooled to very low temperatures. A sharp peak in the velocity distribution was observed 
below a critical temperature, indicating that condensation has occurred, with the alkali 
atoms condensed in the same ground state and showing a narrow peak in the momentum 
space and in the coordinate space J4]. The Bose-Einstein condensation of hydrogen atoms 
|5 1, of calcium atoms js], and of strontium atoms ?! have also been observed. 

Cosmological data provide compelling evidence that about 95% of the content of the 
Universe resides in two unknown forms of energy that we call dark matter and dark energy: 
the first residing in bound objects as non-luminous matter js], the latter in the form of a 
zero-point energy that pervades the whole Universe The dark matter is thought to be 
composed of cold neutral weakly interacting massive particles, beyond those existing in the 
Standard Model of Particle Physics, and not yet detected in accelerators or in dedicated 
direct and indirect searches. There are many possible candidates for dark matter, the most 
popular ones being the axions and the weakly interacting massive particles (WIMP) (for a 



review of the particle physics aspects of dark matter see 1Q|). Their interaction cross section 



with normal baryonic matter, while extremely small, are expected to be non-zero and we 
may expect to detect them directly. Scalar fields or other long range coherent fields coupled 



to gravity have also intensively been used to model galactic dark matter 



111- 



At galactic scales ~ 10 kpc, the ACDM model meets with severe difficulties in explaining 
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the observed distribution of the invisible matter around the luminous one. In fact, N- 
body simulations, performed in this scenario, predict that bound halos surrounding galaxies 
must have very characteristic density profiles that feature a well pronounced central cusp, 
PNFwir) — Ps/( r / r s)(l + r / r s) 2 ll2|, where r s is a scale radius and p s is a characteristic 
density. On the observational side, high-resolution rotation curves show, instead, that the 
actual distribution of dark matter is much shallower than the above, and it presents a nearly 
constant density core: ps(r) = po r o/( r + r o)( r2 + r o) |13(] , where r is the core radius and 
Po is the central density. 

The possibility that the dark matter, which is required to explain the dynamics of the 
neutral hydrogen clouds at large distances from the galactic center, could be in the form 
of a Bose-Einstein condensate was proposed initially in To study the condensate the 
non-relativistic Gross-Pitaevskii equation was investigated numerically. An alternative ap- 
proach was introduced in jl5|. By introducing the Madelung representation of the wave 
function, the dynamics of the system can be formulated in terms of the continuity equa- 
tion and of the hydrodynamic Euler equations. Hence dark matter can be described as a 
non-relativistic, Newtonian Bose-Einstein gravitational condensate gas, whose density and 
pressure are related by a barotropic equation of state. In the case of a condensate with quar- 
tic non-linearity, the equation of state is polytropic with index n = 1. To test the validity of 
the model the Newtonian tangential velocity equation of the model was fitted with a sample 
of rotation curves of low surface brightness and dwarf galaxies, respectively. A very good 
agreement was found between the theoretical rotation curves and the observational data for 
the low surface brightness galaxies. Therefore dark matter halos can be described as an 
assembly of light individual bosons that acquire a repulsive interaction by occupying the 
same ground energy state. That prevents gravity from forming the central density cusps. 
The condensate particle is light enough to naturally form condensates of very small masses 
that later may coalesce, forming the structures of the Universe in a similar way than the 
hierarchical clustering of the bottom-up CDM picture. Then, at large scales, BEC perfectly 
mimic an ensemble of cold particles, while at small scales quantum mechanics drives the 
mass distribution. Different cosmological and astrophysical properties of condensed dark 
matter halos have been extensively investigated 16l4l8|. 

All of the previous researches have been done by assuming that the condensed dark 
matter is at zero temperature. This assumptions is certainly a very good approximation for 
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the description of dark matter in thermodynamic equilibrium with the cosmic microwave 
background, and for the analysis of the galactic rotation curves. It is already well established 
in condensed matter physics that the zero-temperature Gross-Pitaevskii equation gives an 
excellent quantitative descriptions of the Bose-Einstein condensates for T < O.BTbec, where 
Tbec is the Bose-Einstein transition temperature. This condition is obviously satisfied by the 
dark matter halos of the low-redshift galaxies. However, in the early Universe, immediately 
after the condensation, finite temperature effects could have played an important role, and 
significantly affect the cosmological evolution. 

It is the purpose of the present paper to study the finite temperature properties of the 
gravitationally self-bound Bose-Einstein dark matter condensates of collisionless particles. 
In particular, we focus on the description and the static properties of the condensates in- 
teracting with a thermal cloud. The basic equations describing the gravitationally trapped 
condensate at finite temperature are obtained, and their static solution, describing a con- 
densate in thermal equilibrium with a thermal cloud, is obtained by using the Hartree-Fock- 
Bogoliubov and Thomas-Fermi approximations, respectively. The temperature dependent 
dark matter density and mass profiles are obtained, as well the (temperature dependent) 
radius of the condensate and thermal cloud system. By using the explicit form of the phys- 
ical parameters it immediately follows that the finite temperature effects do not play any 
significant role in the characterization of the physical properties of the dark matter halos at 
low redshifts, for which the static solutions of the standard Gross-Pitaevskii equation jl, 2| 
give an excellent description. 

The present paper is organized as follows. The basic equations describing finite tem- 
perature Bose-Einstein condensates are presented in Section [III The effects of the finite 
temperature on the dark matter halo profiles, trapped by their gravitational potential, are 
presented in Section IIHI We discuss and conclude our results in Section (TV] 



II. BOSE-EINSTEIN CONDENSATES AT FINITE TEMPERATURE 



The Heisenberg equation of motion for the quantum field operator $ descri 



dynamics of a Bose-Einstein condensate at arbitrary temperatures is given by 19 



Ding the 



2ii 



(t,r) 

lh —dT- 



h 2 

-—A + mV grav (f) + (t, r) $ (t, r) 
2m 



$ (t, r) 



(1) 
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where m is the mass of the condensed particle, V grav (f) is the gravitational trapping poten- 
tial, and g' = 4nah 2 /m, with a the s-wave scattering length. Eq. ([1]) is obtained under the as- 
sumption that the interaction potential can be represented as a zero-range pseudo-potential 
of strength g'. The equation for the condensate wave-function \& (t, r) = (&(t,r)\ is ob- 
tained by taking the average of Eq. ([T]) with respect to a broken symmetry non-equilibrium 
ensemble in which the quantum field operator takes a non-zero expectation value. Then for 
the exact equation of motion of \& (t, r) we find 



(t,r) 



h 2 

--A + mVgrav (r) 

2m 



V (t, r) + g' ( $ + (t, f) $ (t, r) $ (t, f] 



(2) 



By introducing the non-condensate field operator ijj (t,r) we separate out the condensate 
component of the quantum field operator, so that [19M21 ] 



$ (t, r) = * (t, r)+ip (t, r) 



(3) 

where the average value of if) (t, r) is zero, ( if) (t, r) ) = 0. By taking into account the 
relations 



$ + (t, f) $ (t, f) $ (t, r) = |*| * + 2 |^| 2 if) + * 2 ^ + + + 2*^ + ^ + ^ + ^, (4) 



we obtain 



$+ (t, f) $ (t, r) $ (t, f) > = — Pc m + + p™** 



ifj + if)if) 



(5) 



where we have introduced the local condensate mass density p c (t, f) = mn c (t, f) = 
rn |* (t,r)\ , the non-condensate mass density p(t,r) = mn(t,r) = m (ip + (t,r) if) (t,r)), 
and the off-diagonal (anomalous) mass density p^h = mm(t,r) = m (t,r)if> (t,r)\ , re- 
spectively. With the use of Eq. (jSJ) in Eq. (j2J) we obtain the equation of motion for * as 
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2l| 



3* (t,r) 
dt 



ft 2 

-—A + mVg mt , (f) + gp c (t, r) + 2gp (t, r) 
2m 



gp^* + g' U + (t, r) if) (t, r) if) (t, r, 



^ (t, r) 



(6) 



where we have denoted g = 47iah 2 /m 2 . Besides the condensate density p c (t,f) the gen- 
eralized Gross-Pitaevskii Eq. fl6]) also contains the non-condensate density p(t,r), the off- 
diagonal non-condensate density pm, and the three-field correlation function yif> + if)if>^, both 
of which have nonzero expectation values due to the assumed Bose broken symmetry. 



The generalized Gross-Pitaevski equation can be transformed to a hydrodynamic 
form by introducing the Madelung representation of the wave function as \l/ (t, f) = 



/p^exp \{i/h) S (t, r)]. Then it follows that Eq. ([6]) is equivalent to the system [19 
dp 



2l| 



dt 



+ V • (p c v c 



2|lm 
n 

dS_ 
~dt ~ 



1 



r 



12- 



Pc 



-mvl 
2 c 



(7) 



(8) 



where the local velocity of the condensate is given by v c (t,r) = (li/m) VS and we have 
introduced a new function r 12 , which plays the role of the local source term in the condensate 
continuity equation. The chemical potential of the condensate is defined as 

h 2 A^r c 



Pc 



2m ^fp, 

where the dissipation term R(t, r) is given by 



+ mV grav (f) + gp c (t, f) + 2gp (t, f) + ReR (t, r) 



(9) 



R(t,r) = — [(**) 2 p™. + (^ + # 
pc L 



(10) 



Eq. 



can be reformulated as the Euler equation of fluid dynamics for the condensate, 
dv r 



ni- 



di 



m 



dv c ._ 
— + {v c ■ V) v c 



-V/i c 



(11) 



There are a number of approximations that can be introduced in Eq. ([6]), and which 
have been investigated in detail in the literature. If all the particles are in the condensate, 
p(t,f) = 0, and the anomalous correlations/?^ and ^ + ^-^^) are absent, Eq. reduces to 
the well-known Gross-Pitaevskii equation [2|. 

The Gross-Pitaevskii equation gives a very good description of the condensate if T << 
Tbec- In fact the Gross-Pitaevskii equation gives a good quantitative description of the 
condensate for T -< 0.5Tbec, a temperature range for which the interaction induced de- 
pletion of the condensate is estimated to be of the order of a few percent 2l|, |22|. The 
physical reason for this possibility is that in the limit where the modes of the system are 
highly occupied (iV fc >> 1), the classical fluctuations of the field overwhelm the quantum 
fluctuations. Therefore these modes may be represented by a coherent wave function. This 
is analogous to the situation in laser physics, where the highly occupied laser modes can be 
well described by classical equations. 
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In the Hartree-Fock-Bogoliubov (HFB) approximation one ignores the three-field corre- 
lation function and the p c (t,r), p(t,r), and p m fluctuations are calculated self- 
consistently 23[. The dynamic Popov approximation consists in ignoring both (?p + ip?p 



and p^ 24j. In the static Popov approximation the fluctuations of the density of the ther- 



mal cloud are ignored by assuming that the non-condensate is always in static thermal 
equilibrium, so that p (t, r) ~ p(f) 

In the static case the anomalous density pm can be obtained as p m = 
—gp c 0-/V) J [(l + 2/p) /2Ep\ dp/ (2nh) 3 , where E p is the local excitation energy (per unit 
mass), V is the volume of the system, and /° is the equilibrium Bose-Einstein distribution 



function 



191 ] . The anomalous density, as well as the three-field correlation function, also de- 



pends on the volume of the system. Since the integral takes values of the order of unity, one 
can approximate the anomalous density as p m ~ / (kT/E p ) gp c / (VE P ), where / (kT/E p ) is 
a function of the temperature with numerical values of the order of unity or less. Therefore 
the contribution of the anomalous density term in the generalized Gross-Pitaevskii equa- 
tion Eq. ([6]) is of second order in the coupling constant g < 1, and its contribution to 
the total density can generally be neglected with respect to the condensate density p c . In 
the limit of high temperatures p m / p c ~ (g/VE p ) (kT/E p ), while in the limit of T — > 0, 
Pm/ Pc ~ {9/VEp). In a low temperature system and in laboratory conditions the contribu- 
tion of the anomalous density to the total density is of the order of a few percents. However, 
once T approaches the critical temperature T c , the anomalous density can play a significant 
role in the dynamics of the condensate. The three-field correlation function has the same 
order of magnitude as the anomalous density, thus also giving a second order correction in 
the generalized Gross-Pitaevskii equation. In the limit of V — > 00, a situation that represent 
an excellent approximation in the case of dark matter halos, with radii of the order of a few 
tenth of kiloparsecs, both the anomalous density and the three-field correlation function 
tend rigourously to zero. 



III. FINITE TEMPERATURE BOSE-EINSTEIN CONDENSED DARK MATTER 
DENSITY PROFILES 

In the following we will restrict our analysis to the range of finite temperatures where the 
dominant thermal excitations can be approximated as high energy non-condensed particles 
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moving in a self-consistent Hartree-Fock mean field, with local energy 19M21| 



e P (t, r) = 7^ + mV g rav (0 + 2g [p c (t, f) + p (t, r}] = + U eff (t, r) , (12) 

where U e ff (t, r) = mV grav (r) + 2g [p c (t, r) + p (t, r)]. Hence in the present approximation 
we neglect the effects of the mean field associated with the anomalous density p^ and the 
three-field correlation function (^p + ijj ; ip'^ . 

Moreover, we consider the case in which the thermal cloud and the condensate are in 
static equilibrium. Therefore there are no stationary equilibrium currents, v c = 0, and the 
source and dissipation terms in the hydrodynamic equations of motion vanish, T 12 = 0, 
R = 0. 

In order to describe the gravitational properties of the Bose-Einstein condensed dark mat- 
ter halo we consider the gravitational interaction within a mean field approximation, by in- 
troducing the self-consistent gravitational potential as V grav (t,r) — J G (r — r r ) |\P (t, r) | 2 df, 
where G (r — f) = — Gm 2 /\f — r"| is the usual gravitational interaction potential. There- 
fore the Gross-Pitaevskii equation describes a self-gravitating Bose-Einstein condensate 
with short range interactions. The gravitational potential V grav (t, r) is determined self- 
consistently by the Poisson equation, given by 

AV grav = 47rG[ Pc (f)+p(r)}. (13) 

Under these assumption, and by representing the equilibrium condensate wave function 
as ^ (t,f) = (r) exp [—(i/h)p c t], p c =constant, it follows that the stationary static con- 
densate wave function \l>o (f) satisfies the equation 

h 2 

-—A + mV grav (f) + gp c (f) + 2gp (f) # (r) = (r) ■ (14) 

In the static case the energy of the thermal excitations is given by 

e P (r) = ^ + mV grav (f) + 2g [p c (f) + p (r)} = |- + U eff (f) , (15) 

where the gravitational potential is modified by the Hartree-Fock mean field, U e ff (r) = 
mV grav (f) + 2g [p c (f) + p (f)). 

The collision between particles in the thermal cloud forces a non-equilibrium distribution 
to evolve to the static Bose-Einstein distribution f° (r, p). Hence the particles in the ther- 
mal cloud are in thermodynamic equilibrium among themselves. By using a single-particle 

8 



representation spectrum the equilibrium distribution of the thermal cloud is given by 

/°(f,£)= [ e ^-A_ 1 ]- 1 ) (16) 

where /3 = 1/ksT, with ks Boltzmann's constant, and jl is the chemical potential of the 
thermal cloud. In order to determine ji we assume that the condensate and the thermal cloud 
components are in local diffusive equilibrium with respect to each other. The requirement 
of a static diffusive equilibrium between the cloud and the condensate imposes the condition 
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2l| 



H c = jl. (17) 

Therefore the chemical potential of the condensate also determines the static equilibrium 
distribution of the particles in the cloud. 

As a last step in our analysis we introduce the Thomas-Fermi approximation for the 
condensate wave function. In the Thomas-Fermi approximation, the kinetic energy term 
— (h 2 /2m) A of the condensate particles is neglected. In order to analyze the validity of 
the Thomas-Fermi approximation we consider a system of N bosons, all in the same state, 
occupying a volume V, with spatial extent R, and confined by a gravitational potential. The 
mass of the system is M, while for the mass and the scattering length of the dark matter 
particles we adopt the values m = 3.19 x 10 -37 g and a = 1.78 x 10~ 19 cm, respectively (see 
Section [TV}. The total energy E of the system can be written as E = E^ in + E int + E grav , 
where Ek in , E int and E grav are the kinetic energy, interaction energy, and gravitational 
energy, respectively. The kinetic energy per particle is h 2 /2mR 2 |2|, and therefore the 
total kinetic energy of the system is given by E^ n = Nh 2 /2mR 2 . The interaction energy 
can be obtained as E int = (1/2) (N 2 /V) mg |2|, while the gravitational potential energy is 
Eg^v = GM 2 /R. Therefore the total energy of the system can be written as 

n Ar h 2 3 Ar2 h 2 a GM 2 . . 

E = N ^ + -iV 2 -+ . 18 

2mR 2 2 mi? 3 R V ; 

The interaction energy is much larger than the kinetic energy, E int » E^ in if the condi- 
tion 

3Na 

— » 1 (19) 

is satisfied. By assuming that the mass of the condensate is M = 10 10 M Q = 2 x 10 43 g, 
the dark matter particle number is of the order N = M/m = 6.2 x 10 79 particles. For a 
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condensate with radius R = 10 kpc = 3 x 10 22 cm, the quantity 3Na/R = 1.11 x 10 39 , which 
is obviously much greater than one, 3Na/R » 1. Therefore the contribution of the kinetic 
energy of the particles can be neglected with respect to the interaction energy of the bosons. 
Hence in the case of a large particle number, the Thomas-Fermi approximation gives an 
excellent description of the condensate, at such a precision level that the present description 
of the dark matter halos can be considered exact. The Thomas-Fermi approximation breaks 
down at length scales of the order of R ~ ^/m/47rpa, where p is the mean dark matter 
density. For p = 10" 24 g/cm 3 , we obtain R « 377.65 cm, a length scale that is irrelevant 
from an astronomical point of view. The gravitational energy of the dark matter halo has 
the value E grav = 8.89 x 10 56 ergs, why for the considered parameters E int = 1.35 x 10 56 ergs, 
which is of the same order of magnitude as the gravitational energy. However, E grav > E int , 
thus showing that the gravitational energy is indeed a trapping energy for the dark matter 
halo. 

Hence we can neglect, with a very good approximation, the kinetic energy term in 
Eq. (jUj) . and thus the chemical potential of the condensate (and of the thermal cloud) 
is given by 

p c = p = mV grav (f) + g [p c (f) + 2p (f)] . (20) 
The equilibrium non-condensate density is obtained by integrating the equilibrium Bose 



Einstein distribution over the momentum. Thus we obtain 
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2l| 



777 / 777 

P (r) = ^—3 J d 3 Pf° (r, P) = ^3/2 [z (01 , (21) 

where Xt = ^/2irh 2 (3 /m is the de Broglie thermal wavelength, gs/2(z) is a Bose-Einstein 
function, and the fugacity z (r) is given by 

j[p,-U eff (r)] _ e -Pg Pc {r)_ (22) 



z[r) = e r 



By applying the operator A on both sides of Eq. (|20|) . representing the conservation of 
the energy, and by taking into account that p c and p are constants, we obtain immediately 

gA [p c (f) + 2p (f)] = -mAV grav = -ArrGm [p c (r) + p (f)] . (23) 

Therefore the static condensate density profile at non-zero temperature is described by the 
following equation 

A {p c (f) + 2^,3/2 [e-f*V] } = { Pc (0 + |^3/2 [e-"»>W] } . (24) 

10 



In order to give an approximate representation of the Bose-Einstein functions 
#3/2 {exp [— f3gp c (f)]} appearing in Eq. f l24|) we estimate first the numerical values of the 
quantity —(3gp c (f). As a first numerical approximation we take for the density profile its 
maximum value at the center, so that p c (r) « p c (0). For realistic galactic dark matter 



na. 



os the central density of the dark matter is in the range of p c (0) ~ 10~ 24 — 10~ 26 g/cm 3 
lOt ] . By adopting for the mass and for the scattering length of the dark matter particles 
the values m = 3.19 x 10~ 37 g and a = 1.78 x 10~ 19 cm, respectively, we obtain an upper 
bound for (3gp c (0) as {m&xf3gp c (0)} G (1.7695 x 1(T 7 /T, 1.7695 x 1(T 9 /T). Hence in any 
physically realistic dark matter scenario in which the temperature of the dark matter halo is 
higher than the temperature of the Cosmic Microwave Background Radiation, T > 2.73 K, 
the quantity f3gp c (r) << 1, for all f* and T. As a next step we obtain a linear approxima- 
tion of the function g^/2 {exp [— figp c {f)]}- In order to estimate g3/2{x) we make two linear 
approximations simultaneously, namely, 

9s/2(x) « 53/2(1) + g 1/2 (1) x Ax, (25) 

with x = exp [— f3gp c (f)], and 

Ax = e-to»tT) _ 1 ~ _p gpe (f) . (26) 

Therefore in the first order of approximation the Bose-Einstein function 
<7 3 /2 {exp [— f3gp c (r)]} can be obtained as 

where C (3/2) = (73/2(1) = 2.612 and £ (1/2) = #1/2(1) = —1.460 are Riemann zeta functions. 
Hence the density profile of the condensate component of the finite temperature dark matter 
halo can be obtained by solving the second order differential equation given by 

AirGm 1 - C (1/2) mgf3/X^ 4ttC(3 /2)Gm 2 

Pc[r> ~ 9 l-2C(l/2)m^/A? r Parj ^[l^Ca^m^/A 3 ]' 1 } 

By representing the condensate density as 

Pc (r) = p?> + p? (0 , (29) 

where 

„(o) = C(3/2)m 

Pc A 3 T [l-C(l/2)m^/A|,]' lJUj 
11 



it follows that the function pi (r) satisfies the equation 



ApW (f) = -K* P W 



(31) 



where we have denoted 



K 



Gm 3 1 - C (1/2) mgf3/\ 3 T 



ah 2 I -2C, (1/2) mg(3/\ 3 T 
In static spherical symmetry the general non-singular solution of Eq. ( 131]) is given by 

sin Kr 



(32) 



P « (r) = a: 



(33) 



where A is an arbitrary integration constant. The density of the non-condensate component 
of the dark matter can be written as 



1\ mgf3 



while the total density of the finite temperature dark matter halo is given by 



Ptot{r) = p c {r) +p(r) 



C 



1\ mg(3 



2 J Aji 



sin Kr 
Kr 



(34) 



(35) 



At the center of the halo the density of the dark matter is ptot(fy = which allows us 
to determine the integration constant A as 



.4 



2 / A?r 



(o) 

Plot 



C(3/2)m 



A|[l-2C(l/2)m^/9/A|] A 



m / 3 



The total density profile of the dark matter halo can be represented in the form 



( o)SinKr 

Ptot (r) = Plot— [? — + 
Kr 



m 



P 



(o) 



sin Kr 
Kr 



(36) 



(37) 



The previous results can be written in a more transparent form if we introduce the 
condensation temperature T B ec, given by T BEC = 2irh 2 pf EC /( 2 / 3 (3/2) m 5 / 3 k Bl where Pbec 
is the density of the dark matter at the condensation moment. Then we obtain immediately 



m 

A^T 



rp \ 3/2 

PBEC 1 T * ^/3 = 2C 2/3 (3/2) 



and 



C(3/2) \ Tbec) 
mg(3 



t y 1 



1/3 



Aji 



^(3/2) w>l 3pBEC \T, 



T 



EEC 



1/2 



respectively. Hence we obtain 



K 2 (T) 



Gm 3 l-a (a/m 1 / 3 ) p% c (T/T BEC ) 1/2 
ah 2 l-2a(a/wM 3 )pf EC (T/T BEC ) 1 ' 2 ' 



(38) 



(39) 



(40) 
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where we have denoted a = 2C(l/2)/C 1/3 (3/2) < 0. By using a series expansion for 
( (1/2) mgft/Xj, « 1, K 2 (T) can be approximated as 



K 2 {T) 



Gm 3 
ah 2 



/3 



or 



if 2 (T) 



Gm 3 
a/?. 2 



1i a 1/3 

1 + a— tttPbbc 



1/3' 



T 



1/2' 



(41) 



(42) 



For the total density profile of the dark matter halo we obtain 



Ptot{r) 



(o) 

Ptot 



a 



a (a/m 1 / 3 ) pj EC (T/Tbec) 
l-a(a/m 1 / 3 )pT EC (T/T BEC ) 112 
(a/m 1 / 3 ) p f Ec (t/Tbec) 2 



sin [K(T)r] 
K(T)r 



or, equivalent ly, 



I- a(a/m 1 / 3 )pf EC (T/Tbec) 1 / 2 ' 

sin [K(T)r] 



Ptotir) 



(o) 

Ptoi 



Pt(T) 



where we have denoted 

Pt(T) 



a 



K(T)r 
(a/m 1 / 3 ) p f Ec (t/Tbec) 



-Pt(T), 



< 0. 



l-a(a/m^)pf EC (T/T B Ec) 1/2 
The mass profile M to t(r) = An f Q r ptot(r)r 2 dr of the dark matter halo is given by 

Ana (a/m 1 / 3 ) pf EC (T/T BEC f 



(43) 



(44) 



(45) 



M tot (r) 



a(a/mV3 )p ^ c (T/T r - 



K 3 (T) 



(o) 

Ptot 



a 



BSC J 
,V3^| n 4 Jl. 



■r 3 + 



(a/mV3) p % EC {T/Tbec) 



a (a/m 1 / 3 ) p^ EC (T/Tbec) 1 



/2 



x 



{sin [AT(T)r] - K(T)r cos [X(T)r]} 



(46) 



or, equivalently, 
An 



Mi 



tot 



3 Pr(T)r3 + #^t) 



p T {T) {sin[K(T)r]-K(T)r cos [K(T)r]}. (47) 



The radius R of the dark matter distribution is determined by the condition ptot{R) — 
In the zero-temperature limit T — > 0, the radius i?o of the halo is given by Ro = it/K(0) 
By representing K(T)R as K(T)R = n + </>ttt, we obtain first sm[K(T)R] = — sin (<^7r) » 
— 0r7r. Hence for 0^ we obtain 

Pt(T) 



MT) 



p!2+2p T (T)■ 



(48) 
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Therefore the radius of the dark matter halo at a finite temperature T is given by 



R(T) 



7T 



K{T) 



Pt{T) 
p£l + 2p T (T) 



(49) 



The tangential velocity of test particles in stable circular orbits in the finite temperature 
dark matter halos can be obtained as v^ g = GM to t/r. The variations with respect to the 
dimensionless radial coordinate 9 = K(T)r of the dark matter halo density and mass are 
presented, for different values of the temperature, in Figs. [1] and [2], respectively, while the 
temperature variation of the radius of the dark halo is shown in Fig. [3l 



o o 
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0.8 



0.6 



0.4 



0.2 



0.0 








4 







FIG. 1: Radial coordinate dependence of the total density profile of the finite temperature dark 
matter halo, for different values of the temperature: T = (solid curve), T/Tbec = 0.2 (dotted 
curve), T/Tbec = 0.4 (dashed curve), T/Tbec = 0.5 (long dashed curve), and T/Tbec = 0.6 
(ultra-long dashed curve). For all the curves a (a/m 1 / 3 ) = —5000 cm/g 1 / 3 . The Bose-Einstein 
transition density is taken as pbec = 10~ 21 g/cm 3 , while the central density of the dark matter 
halo is pfo t = 10~ 24 g/cm 3 . The dimensionless radial coordinate 9 is defined as 9 = K(T)r. 
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FIG. 2: Radial coordinate dependence of the total mass profile of the finite temperature dark 
matter halo, for different values of the temperature: T = (solid curve), T/Tbec = 0.2 (dotted 
curve), T/Tbec = 0.4 (dashed curve), T/Tbec = 0.5 (long dashed curve), and T/Tbec = 0.6 
(ultra-long dashed curve). For all the curves a (a/rra 1 / 3 ) = —5000 cm/g 1 / 3 . The Bose-Einstein 
transition density is taken as pbec = 10 -21 g/cm 3 , while the central density of the dark matter 
halo is pfo t = 10~ 24 g/cm 3 . The dimensionless radial coordinate 9 is defined as 9 = K(T)r. 

IV. DISCUSSIONS AND FINAL REMARKS 

In the present paper we have analyzed the effects of the finite temperature on the Bose- 
Einstein dark matter halos. We have obtained the equations describing the condensate and 
the non-condensate component of the dark matter at finite temperatures T > 0, and we 
have derived the density profile, the mass distribution and the radius of the static finite 
temperature configuration in thermodynamic equilibrium. To obtain the description of the 
finite temperature dark matter halos we have considered an analytical approach, in which 
we have ignored the anomalous density term, and used a phenomenological description for 
the thermal excitations and non-condensed particles. This approach, despite its simplicity, 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 

T/7bec 



FIG. 3: Radius-temperature relation for the finite temperature dark matter halos. The values of 
the physical parameters are a (a/m 1 / 3 ) = —5000 cm/g 1 / 3 , pbec = 10~ 21 g/cm 3 , while the central 
density of the dark matter halo is pfj t = 10~ 24 g/cm 3 . 

gives a very clear physical description of the finite temperature properties of the dark matter 
halos. 

We have explicitly determined the temperature dependence of the density profile of the 
dark matter halo, given by Eq. As one can immediately see, if T «Tbec, t ne zero 

temperature density profile can give an excellent description of the condensate in thermo- 
dynamic equilibrium with a thermal cloud. The effect of the presence of the thermal cloud 



results in the modification of the central density of the system [25| , which consists now from 
the sum of the central density of the condensate and of the thermal cloud. However, in order 
to see if the approximation of the zero temperature indeed works, one have to estimate the 
temperature at which the Bose-Einstein condensation took place, and the corresponding 
density of the Universe. 

We assume that in the early stages of the evolution of the Universe dark matter consisted 
of bosonic particles of mass m and temperature T, originating in equilibrium, and decoupling 
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at a temperature To or chemical potential [ir> » m. By assuming that the dark matter 
forms an isotropic gas of particles in kinetic equilibrium, the pressure of the non-condensed 
gas can be obtained as [17]] 

P = PDMC 2 a 2 v , (50) 

where pdm is the density of the dark matter particles before condensation, a v = (w 2 )/3c 2 , 
and (v 2 ) is the average squared velocity of the particles, ay is the one-dimensional velocity 
dispersion. By taking into account the continuity of the pressure during the condensation 
phase one can determine the transition density Pbec- The numerical values of the physical 
parameters at the transition point depends on three unknown parameters, the dark matter 
particle mass m, the scattering length a, and the dark matter particles velocity dispersion 
ay, respectively. By assuming a typical mass of the dark matter particle of the order of 1 
eV (1 eV = 1.78 x 10~ 33 g), a typical scattering length of the order of 10~ 10 cm, and a mean 
velocity square of the order of (v 2 ) = 81 x 10 14 cm 2 /s 2 , the critical transition density can be 



written as 



p BEC = 3.868 x 10 



-21 



ar 



3 x 10- 6 



X 



m 



10-33 gj VlO" 10 



cm 



g/cm 3 . (51) 



by 



he critical temperature at the moment of Bose-Einstein condensate transition is given 



Tbec ~ 6.57089 x 10 3 x 



m 



10- 33 g 



1/3 



a 



V 



3 x 10" 



2/3 



10 



-10 



cm 



-2/3 



K. 



(52) 



For the critical redshift we obtain the value 171 ] 



l + z c 



1.224 x 10 3 x 



m 



<7t 



3 x 10- 



10 -33 
l/3(l+4) 



l+a( 



X 



10 



-10 



cm 



-l/3(l+(T 



(53) 



In obtaining these values for the critical for Pbec-, Tbec an d z cr we have implicitly 
assumed that the Bose-Einstein condensation process took place relatively late in the history 
of the Universe, at a redshift z cr w 1000. For a discussion of the possibilities of a Bose- 
Einstein condensation in the very early Universe see 26|. By assuming that dark matter 
is cold, and in thermal equilibrium with the cosmic microwave background, it follows that 
the thermal contribution to the central density of the condensate is at least four orders 
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of magnitude smaller than the central density of the dark matter, which may be of the 
order of 10~ 24 g/cm 3 . Hence these simple qualitative estimation also show that for low 
redshift galaxies the general form of the density profile of the Bose-Einstein condensate dark 
matter can be described within an excellent approximation by the zero-temperature profile, 
p = p c sin(Kr)/Kr, even if we assume that the Bose-Einstein condensation of dark matter 
took place relatively late during the cosmological evolution of the Universe. 

The zero-temperature approximation also allows us to make an estimate of the dark 
matter particle mass. At the boundary of the dark matter distribution p(R) = 0, giving 
the condition KR = n, which fixes the radius of the condensate dark matter halo as R = 
7r^h 2 a/Gm 3 . The total mass of the condensate dark matter halo M can be obtained as 
M = 4vr 2 (h 2 a/Gm 3 f /2 p c = AR?p c /it, giving for the mean value < p > of the condensate 
density the expression < p >= 3p c /ii 2 . The mass of the particle in the condensate is given 
by [15| 

(2 fc2 \ 1/3 
J « 6.73 x 1(T 2 x [a (fm)] 1/3 [R (kpc)]^ 2/3 eV. (54) 

For a ~ 1 fm and R ps 10 kpc, the typical mass of the condensate particle is of the order 
of m ~ 14 meV. For a ~ 10 6 fm, corresponding to the values of a observed in terrestrial 
laboratory experiments, m ~ 1.44 eV. These values are consistent with the limit m < 1.87 
eV obtained for the mass of the condensate particle from cosmological considerations 26 ]. 

The properties of dark matter can be obtained observationally from the study of the 
collisions between clusters of galaxies, like the bullet cluster (IE 0657-56) and the baby 
bullet (MACSJ0025-12). From these studies one can obtain constraints on the properties of 
dark matter, such as its interaction cross-section with baryonic matter and the dark matter- 
dark matter self-interaction cross section. If the ratio a m = a/m of the self- interaction cross 
section a = 4ira 2 and of the dark matter particle mass m is known from observations, then 
the mass of the dark matter particle in the Bose-Einstein condensate can be obtained as 



,^h 2 ^\ 2 " 
m = jyT ■ ( 55 ) 



2G R? t 

By comparing results from X-ray, strong lensing, weak lensing, and optical observations 
with numerical simulations of the merging galaxy cluster IE 0657-56 (the Bullet cluster), 
an_upper limit (68 % confidence) for a m of the order of a m < 1.25 cm 2 /g was obtained in 
. By adopting for o m a value of a m = 1.25 cm 2 /g, we obtain for the mass of the dark 
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matter particle an upper limit of the order 



m < 3.1933 x 10 



-H7 



R V 4/ V a m 



10 kpcy \1.25cm 2 /g 
R V 4/5 / a„ ^ 1/5 



0.1791 x — ^— meV. (56) 

V 10 kpcy \1.25cm 2 /g/ 

By using this value of the particle mass we can estimate the scattering length a as 



a < \l am ^J n = L7827 x 10 ^ 19 cm = L7827 x 10 ^ 6 fm - ( 5T ) 

This value of the scattering length a, obtained from the observations of the Bullet cluster 
IE 0657-56 is much smaller than the value of a = 10 4 — 10 6 fm corresponding to the BEC's 
obtained in laboratory terrestrial experiments 

If Bose-Einstein condensate dark matter is formed of Cold Dark Matter elementary parti- 
cles, their superfluidity properties may distinguish them from other dark matter candidates 
due to the presence of quantum vortices. The presence of vortices in the central part of 
the BEC dark matter halos with their empty core change the gravitational coupling with 
baryons and contribute to a nearly smooth dark matter distribution in the halo. The vortex 
stability condition is satisfied when the rotating halo has enough angular momentum to sus- 
tain more than one vortex like formation of vortex lattice due to rotation of spiral galaxies. 
Because of this rotation vortices manifesting as density holes appear in the condensate with 
quantized circulation, V = § c v ■ df = I (h/m), where T is the circulation, C is any contour 



around a vortex, and I is the topological charge [27J. T is quantized in units of h/m. When 
I = 0, there is no vortex. A rotating environment means that / = 1 with one stable vortex or 
for the case of / > 1 unstable vortices, which decay to the case having 1 = 1. Thus, rotation 
of a self-gravitating dark matter condensate is a source for quantized vortices, which evolve 
towards a vortex lattice. Its number density varies as n v oc r _1 , and suggests a flat velocity 
profile. 



In a recent paper 



29] , the authors claim to have ruled out the Bose-Einstein condensate 



dark matter model. In order to support this conclusion an equation of state for 



;he finite 



temperature condensate is derived. The isothermal equation of state obtained in [29] indicate 
a Bose-Einstein condensed core surrounded by a non-degenerate envelope, with an abrupt 
density drop marking the boundary between the two phases. By analyzing the behavior of 
the galactic rotation curves the authors conclude that since such a behavior is not confirmed 
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observationally, the Bose-Einstein dark matter condensate model can be ruled out. The 
detailed analysis of the finite temperature effects performed in the present paper has shown 
that once the thermal effects are correctly taken into account, such a conclusion is incorrect. 



The main physical process missing in the analysis in [29| is the requirement of a static 
diffusive equilibrium between the cloud and the condensate, which forces the thermal cloud 
to follow the density distribution of the condensate. Moreover, since the temperature of 
the dark matter in low redshift galaxies is much smaller than the transition temperature, 
thermal effects do not play any significant role in the description of the dark matter halo 



density profiles, a result well-established in the condensed matter literature |19l-l2ll]. This 
situation is very similar to the case of other astrophysical systems having a pure quantum 
nature, like, for example, white dwarfs and neutron stars. Quantum statistical systems in 
which the thermal energy is much lower that the transition energy can be described with a 
high level of precision by the zero temperature limit. For example, the maximum mass and 
the physical properties of the white dwarfs and neutron stars (including the Chandrasekhar 
mass limit) are obtained for the T = case, even that the temperature of the stars may be 
of the order of T = 10 6 — 10 8 K. But the T = approximation gives a very good description 
of the star because the thermal energy of the star is much lower than the Fermi energy of the 
electrons or neutrons. The situation for the dark matter halos is similar-if their temperature 
is much lower than the transition temperature, the T = approximation gives an excellent 
description of the system. 

In conclusion, in the present paper we have shown that for temperatures much lower 
than the Bose-Einstein transition temperature, < T < 0.5Tbec the finite temperature 
effects on the density profiles of the condensed dark matter halos can be neglected. The 
zero-temperature density profiles gives an excellent description of the observed properties 
of the rotation curves [15|. However, for T > 0AT BEC , finite temperature effects can play 



an important role in determining both the static and dynamic properties of the condensate. 
Hence finite temperature effects can play an important role in a cosmological context and 
for the analysis of the cosmological dynamics after the condensation. 
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